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Abstract 

We show how the Fourier transform of a shape in any number of dimensions can be simplified 
using Gauss's law and evaluated explicitly for polygons in two dimensions, polyhedra three 
dimensions, etc. We also show how this combination of Fourier and Gauss can be related 
to numerous classical problems in physics and mathematics. Examples include Fraunhofer 
diffraction patterns, Porods law, Hopfs Umlaufsatz, the isoperimetric inequality and Didos 
problem. We also use this approach to provide an alternative derivation of Davis's extension 
of the Motzkin-Schoenberg formula to polygons in the complex plane. 

1 Introduction 

A shape can be defined mathematically in many ways. Here we consider defining it by a function 
that has the value unity inside the shape and zero outside the shape. Only simply connected 
orientable [1] shapes will be considered. The Fourier tranform of a shape defined in this way has 
interesting applications and connections to many areas of physics and mathematics. For example, 
in physical optics the Fourier transform of the shape of an aperture or opening in an opaque screen 
yields the Fraunhofer or far field diffraction pattern that is generated when light passes through the 
aperture [2]. In X-ray scattering, the Fourier transform of a volumetric shape in three dimensions 
reduces in the appropriate limit to Porods law [3]. From a probabilistic point of view the Fourier 
transform of a shape properly normalized can be considered as the characteristic function or moment 
generating function of the shape. Hence the Fourier transform of the shape is intimately related to 
the moments of the shape and hence to the "shape from moments problem" , which has been studied 
recently for polygons by Golub, Milanfar and others [4] [5] and is related to the overall problem of 
pattern recognition. The Fourier transform of an area bounded by a smooth curve or a polygon 
in the plane can be shown to imply Hopfs Umlaufsatz [6], which states that the tangent vector 
(or the normal vector or any linear combination thereof) to a smooth or piecewise smooth closed 
orientable simply connected curve in the plane rotates by ±2ir as it goes completely around the 
curve. The +/— signs correspond to counterclockwise/clockwise rotation. This seemigly obvious 
geometrical fact is rather subtle to prove [6] We also use this formula to derive the planar version 
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Stokes law [7] and to show that a circle maximizes the area for a given perimeter length which is 
known as the isoperimetric inequality [6] [8]. Finally there is the famous problem faced by Diclo, 
Queen of Carthage, which was to determine what shape open curve of a given length encloses the 
maximum area when its endpoints are connected by a straight line with an arbitrary length [9]. 
Here I show how combining the Fourier transform with Gauss's Law provides a central platform 
for relating and solving all these problems. 

As an aside, the Fourier transform is just one of many techniques that are used for shape 
discrimination [10]. This problem is both mathematically interesting, and, given the ubiquity of 
digital data bases, technologically important. Here our specific interest is in the relation of the 
Fourier transform of a shape to its moments and to the problems mentioned above. This is related 
to, but generally different from, the issues concerned with using Fourier transforms explicitly for 
shape discrimination in digital data bases and we will not discuss this aspect further. 

The paper is organized as follows. In the Section 2, Gauss's Law is used to rewrite the Fourier 
transform of a D dimensional volume as the surface integral over the D — 1 surface of that volume. 
This suface integral can then be evaluated exactly for polygons in two dimensions, polyhedra in 
three dimensions, etc. In particular in two dimensions this then provides an explicit relation between 
the vertices of a polygon and its moments. This section also discusses how this result is related 
to the isoperimetric inequality, Didos problem and Hopf's Umlaufsatz. Finally we show how the 
combination of Fourier and Gauss provides an alternative derivation of a famous result of Davis 
concerning the integral of an analytic fuction over a polygon in the complex plane [11] [4]. Section 3 
briefly discusses how the solution for two dimensional polyhedra appears when Fourier transforming 
three dimensional polyhedrons before showing how the surface integral form of the three dimensional 
Fourier transform can be used to derive not only the standard Porods law [3] for X-ray scattering 
from spherical particles but also the extension of this law to anisotropic particles. This extension 
has been discussed recently in a series of papers by Ciccariello and colleagues [12]. The derivation 
presented here is somewhat more direct than that used by Ciccariello and colleagues. 



2 Fourier and Gauss 

For generality we begin in D dimensions but will rapidly particularize to 2 and 3 dimensions. Define 
a shape, in D dimensions, by the function 9v (x) where 9y (x) = 1 for x = (xi,X2, ■ ■ ■ ,^_d) inside 
V and for x outside. This is known in certain circles as an indicator function [13]. Here we take 
Xi for i = 1, 2, • • • , D to be Cartesian coordinates. The (normalized) moments (x Pl ■ ■ ■ x p ^) of the 
shape are then given by 

« ...a*p) = l J d D x 6 V (x) x Pl • • • x p D D = J J d D x x Pl ■ ■ ■ x p D D (1) 

v 

where J v — J 6y indicates integration over the shape, v = J d D x 0y (x) is the finite volume of the 
shape and pi = 0, 1, 2, 3, ....for each i. The symbols V and v are used to distinguish between the 
shape itself and its volume as a numerical value. Also the terms "volume" and "surface" are used 
generically and should be understood to refer to a D dimensional submanifold in D dimensions and 
to its D — 1 dimensional boundary, respectively. 

Note that (x) can be thought of as a probability density for x to be inside the shape. The 

characteristic function or moment generating function 4>\P) is then given by the Fourier transform 
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9 V {x)/v, i.e., 



d D x 9 V (x) 



(2) 



with 0-x — J2f=i Pi x i = Pi x i- (Unless noted otherwise the Einstein summation convention, wherein 
repeated indices are summed over their appropriate range, will be used throughout the paper.) The 
bar on <f> is meant to indicate that this is the normalized charactistic function, i.e., (0) = 1. Below 

we will find it convenient to work with the unnormalizcd function (f> (jPj . 



Using the power series representation of exp 



i(3 ■ x 



and the multinomial theorem it follows that 



d D x 



E 

n=0 



(iPiXiY 



n=0 P1+P2-I \-pD=n 



PlW ■ ■ -PD ] - 



(3) 



This result indicates the simple but direct relation that exists between the Fourier transform of 
a shape and its moments (x^ 1 ■ ■ ■ x^f ) . Note that we have assumed that since v is finite we can 
interchange orders of integration and summation. For all practical cases this is certainly true. 
We see from this result that if we can compute cj> explicitly then the moments (x^ 1 ■ ■ ■ a^f ) 

follow directly from a Taylor expansion of cj> (jlj about (3 = 0. An explicit form for (jlj can be 
found by applying Gauss law to the Fourier integral. Noting that 



a- 



i/3-x 



(4) 



with d = (d/dxi,d/dx 2 , • • • , d/dx^) = (81,82, • • • , do) we get 

J d D xe^ s = J ' d D xd- (j^^A 



J,- J d D -h^jn(s)e^ 
dV 



(5) 



Here dV indicates the surface of V with s = (s\,S2, ■ ■ ■ ,Sd-i) being coordinates on the surface 
dV . R (s) gives the position in D dimensions of the point on the surface labeled by s and 



g (s) = |det [gij (s)]\ = det d 8t R(s) ■ d Sj R(s) 



(6) 



where (s) is the induced metric on dV, h (s) is the local outward normal to the surface dV and 

d 2 s 



indicates the absolute value [14]. Below we use (5) to prove that 



dct 



d Si R{Z) ■ d Sj R{s) 



is the area element of a parallelogram in 2 dimensions with sides defined by the vectors dRi 
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d Si R(s) dsi (no sum on «).and that ^/|dct [Si ■ dj]\ is the volume of the parallelogram formed by 3 

d 3 s 



noncoplanar (nonzero) vectors Si, for i = 1, 2, 3 in three dimensions and so \j dot d Si R (s) ■ d Sj R (s) 

with i, j = 1, 2, 3 is the volume element in three dimensions. 
Combining (3) and (5) yields 



dV 

= E ^^•••^^••'^) (7) 



n=0 pi+P2H Vpm=n 



It will be convenient below to work with the unnormalized form of the characteristic function 

(£) - v4> ($) (8) 

and write the unnormalized moments as 

" (xf ■ • • x v g ) = v K 1 • • • x p D D > = J d D x 9v (x) xf ■ ■ ■ x p D D (9) 



Ml 



Before particularizing to 2 and 3 dimensions note that expanding the left hand side of (7) in powers 



of/3 = 



yields a term proportional to This term must vanish since there is no corresponding 



power of (3 on the right and we get the result that 

d°- 1 s^{i)n{3) = {) (10) 



/■ 

dV 

This can be interpreted as the statement that "the boundary of a boundary is zero" . This boundary 
of a boundary principle has been considered to have interesting implications with respect to the 
origin of physical laws. [15] 

The remainder of the paper considers particular cases in 2 or 3 dimensions where (j) (jPj can be 

computed from the surface integral either exactly or approximately. We now relate these results to 
various "classical" problems in physics and mathematics. 

3 Results in Two Dimensions 

In this section we show how (5) and (7) can be used to derive various two dimensional classical 
results in physics and mathematics. 

3.1 Smooth Curves 

For smooth curves in two dimensions (5) reduces to 

j d 2 xe^ S = A . J ' ds n( S ) e^M (11) 

V dV 
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where V is the region enclosed by the curve dV given by R(s) with the parameter s being the 
length along the curve. With ds defined as the unit of length on the curve, i.e., ds 2 = dR 2 , it 
follows that g (s) = 1. The unit tangent vector to the curve at s is given by 

t(s) = d s R{s) (12) 

which is automatically normalized since ds = i.e., \d s R (s) = 1. (Generally we will use a hat 

"'" to indicate the vector has been normalized to have unit length. Also we will switch back and 
for between vector and component notation, e.g., between writing u\ and (1,0).) 

We take increasing s to correspond to counterclockwise circulation of the curve so that 



where 



h (s) = £ • i (s) = e ■ d s R (s) 



1 

-1 



(13) 
(14) 

(15) 



and the "•" indicates matrix multiplication, i.e., 

Hi (s) = EijUj (s) 

with repeated indices summed over 1,2. Note that e is the D — 2 version of the totally antisymmetric 
tensor, often called the Levi-Civita tensor [16], defined as +1 for even permutations of i, j = 1,2, 
— 1 for odd permutations and if i = j. 

Substituting (13 ) into (11) and integrating by parts yields 

1 



J d 2 xe^ 



if 



ds 



dV 



= - J ds {p-e-R (*)) (p ■ d s R (s)) e^ s 

dV 



(16) 



The nominal 1//3 pole in (11) has been cancelled and we can set f3 to zero for arbitrary nonzero 



(3 = 13/ 



which gives 



v = J ' d 2 x = - J ds (p-e-R (s)) (p ■ d s R (sfj 



V dV 

Letting ft be (1,0) or (0,1) with R(s) = (Ri {s),R 2 (s)) yields 



v = — J ds R 2 (s) d s R\ (s) = J ds Ri (s) d s R2 (s 

dV dV 



(17) 



(18) 



Averaging the two forms gives 



i J ds (R 1 (s) S R 2 (s) - R 2 (s) d s R± (s)) 

dV 

i J ds Ri (s) Sijd s Rj (s) 



(19) 



Both (18) and (19) are standard forumulae for the area v enclosed by the curve R (s) [17]. 
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Isoperimetric Inequality There are numerous proofs of this inequality [6] [8]. One rather cute 
and straightforward proof uses Fourier Series. Since dV is a closed curve (the boundary of a 
boundary is zero) the functions Ri (s) are periodic in the length 



ds (20) 

dV 

and so can be written as Ri (s) = J2n=-oo R*,n ex P [n^iris / Xjwith R i n = R* _ n so that X t (s) is real 
("*" indicates complex conjugation). Substituting into (5) and using some very simple inequalities 
yields one proof of the isoperimetric inequality (see for example, Luthy in [8]). 

The approach we will use use is to extremize v with the constraint that (j) s R (sfj = 1 imposed 
using a Lagrange multiplier A. Minimizing 

v = ij dsRi ( s ) £ ..Q sRj ( s ) + \Uds (d s R t (s)f - L (21) 
av \o / 

with respect to Ri(s) , and A yields 

01) 

= 7Tr~TT = ( s ) — 2Xd s 2 Ri (s) 
oRi (s) 

Sv 

= = -d s R 1 (s) - 2X8 S 2 R 2 (s) (22) 



Li 



Here 5/5Ri (s) indicates functional differentiation [18]. Combining the first and second equations 
in (22) gives 

(d s Ri) + 2Xd 2 (d s R t ) = for i = 1, 2 (23) 
Taking d s R\ (s = 0) = so that <9 s i?2 (s = 0) = 1, the solutions are 

<9 SJ Ri(s) = -sin(s/2A) 

d s R 2 (s) = cos(s/2A) (24) 
where the minus sign was chosen for convenience. Integrating both sides yields 

i?i(s) = 2Acos(s/2A) +Ro,i 

R 2 (s) = 2Asin(s/2A) +i? ,2 (25) 

which is a circle of radius 2A centered at the arbitrary position Ro — (Rq,i, R0.2) ■ Note that the 
radius 2A is indeterminate since substituting these solutions into the third equation above yields an 
identity. This is as it should be since the circle does not have to be a particular size to maximize an 
unspecified value of v. It just has to be a circle. Of course once a particular value of v is given then 
A can be determined. Note that the above calculation shows only that the circle is one shape which 
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extremizes the area. In other words, we have only shown that the circle is a local extremcum, to 
formally complete the proof requires showing that it is the global extremum. 

We will not show it here but isoperimetric inequality is generalizable to arbitrary dimensions 
with the result that the shape that maximizes the "volume" for a given surface "area" is always a 
"ball" [6] [8], i.e., circle for D = 2, sphere for D = 3, etc. 

The isoperimetric inequality can be used to solve Dido's problem which is to find the open curve 
of fixed length L which encloses the maximum area when the end points are connected by a straight 
line of arbitrary length [9]. Combining this closed curve with its reflection relative to the straight 
line yields a closed curve of length 2L. Via the isoperimetric inequality the circle with circumference 
2L encloses the maximim area and hence the original open curve which maximizes the area is a 
semicircle. 

Stokes Law We now use (5) to derive the planar version of Stokes law as used in classical 
electrodynamics [7]. It should come as no surprise that we can relate Gauss's law to Stokes law 
and indeed to Green's theorem as well since they are all different aspects of the general version of 
Stokes theorem represented compactly using differential forms [16]. Consider the integral over an 
area V in the X1X2 plane of the curl of a smooth 3 dimensional vector field F (x). Let indices from 
the beginning of the aphabet a, 6, c, • • • range over 1,2,3 and indices from the middle of the aphabet 
k, ■ ■ ■ range over 1,2. This way F (x) = F a (x) u a and x — x a u a with u a the unit vectors with 
respect to x a . Represent F (x) as a Fourier transform 

Ft (x) = J d 3 m (£) (26) 

with ft = (3 a u a - Again we assume that orders of integration can be freely exchanged. Using index 
notation the a th component of the curl can be written 

(d*F{x)) a =s ahc d b F c (x) (27) 

Here e a bc is the D = 3 totally antisymmetric Levi-Civita tensor defined as +1 for even permutations 
of a,b,c = 1,2,3, —1 for odd permutations and if any two of the indices have the same value. 
Repeated a, b, c, • • • indices are summed over the range 1,2,3. The normal to V, which lies in the 
xix 2 plane, is U3 and so 

J d 2 x i 3 '(ax P{2 v j) = J d 2 xe 3bc d b F c (x v ) (28) 
v v 

where xy = XiUi. But £3b c vanishes unless b and c are restricted to the range 1,2 and so £sbcdbF c (x) = 
SijdiFj (a?) where Ejk are the elements of the 2x2 Levi-Civita matrix e defined above. Putting this 
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together and using (5) and (26) gives 

J d 2 x u 3 - (dx F(x v fj = J d 2 x EijdiFj (xy) 



= J cPpipiSvFj ($) J d 2 x 
v 

= J ' d 2 W t e l0 F (p)^p- J 'dst,(s)> 

dV 

= J dsU (s) J tfpe&^Fi (p) 

dV 

= j ds i(s)-F^R(s^j 

dV 

which is the planar form of Stokes law. To get the fourth line we used 



iP-R(s) 



[3 1 j ds h (s) , 

dV 



-j3 2 J ds t 2 (s) . 

dV 



which follows from 



ds i/3 -i (s) e 



ifj-R(s) 



dV 



J ds 0-d s R(s)e^-^ 

dV 

J ds d s e^^ 



dV 




3.2 Polygons 



(29) 



(30) 



(31) 



A polygon with N sides and N vertices in two dimensions can be defined by its vertices, arranged in 
a particular order, v\, V2, • • • , vm with Vi = (xu, X2i) ■ For computational convenience it is useful to 
let vn+1 = v±. We consider only orientable polygons, i.e., none of the edges cross over or intersect 
one another. The integral over dV in (7) can be evaluated explicitly in this case with the result 



0) = J*"/ d D -hn(s).e^ 

dV 



R(s) 



N 7? 



where 



-i t % ?v ( cxp ' ~ cxp ' ^)) 

P „=i 0-K+l-«n) V V / V // 

^^^• £ = (/3 1 ,/3 2 ). K Jj =(/3 2 ,-/3 1 ) 



(32) 



(33) 
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The "•" stands for standard matrix vector multiplication. 

Equating the surface integral to the series expansion in terms of moments gives 



E l - E u n ' vM *~ p)M ^ 



n=0 p=0 



£ r£?-7; (- ■ «•») - exp 0) <34) 



£ „ =1 /3 • (Vn+l - V n ) 

Now expand the right hand side in powers of f3 i 



1 ^2 £± K+ 1 ^ CX p (0 . y n+ A _ CXp ^ . tfj ) 

P n~i P-(v n +i-v n ) V V 7 \ // 



m! 

. m— 1 



„=l /? ' («n+l - #n) 

32 E E ~|^-L ' (^+1 _ ^) E {fi ' ' ( 35 ) 



^ m— 1 n— 1 p— 

where 1 — (3 = i i — 0i + /3 2 2 . In the second step we have used the identity 



a n - b n = {a-b)J2 a n - 1 - m b m (36) 



m=0 



Substituting (35) into (34) gives 



oo n 



E l - E ,, ^r p M (x^r p ) 

^ ^!^ P!(«-P)! 

. oo AT m — 1 -i 

= ~~p? E E ' (Vn+l - Vn) 22 (rf -Vn+l) [10 ' V n J (37) 

^ m—1 n—1 p— 

Letting (3 = (3$ and /3 ± = (3$ ± we see that the n th term on the left hand side is proportional to 
f3 n whereas the m th on the right is proportional to f3 m ~ 2 . Hence considering each side as a power 
series in [3 the coefficients of the m term on the right must equal the coefficient of n = m — 2 term 
on the left. Below we write out the general relation between M (xfx^) and powers of the vertices. 
Here we begin by considering the first few terms individually. 

3.2.1 to = 1 term 

The to = 1 term on the right has no corresponding term on the left of (34) and so we must have 

N 

= ' E - *») ( 38 ) 

P n=l 
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The fact that this should vanish follows from the fact that it scales as 1//3 whereas the left hand 
side remains finite as /3 — > 0. Also this term is imaginary and the result must be real. And indeed 
this term does vanish, trivially, since we have defined vn+i = v\ and the sum of the directed sides, 
v n +i — v n , of a closed polygon must vanish. For the case where the boundary dV is everywhere 
smooth the form of the dominant term as (3 — > is given by 

J %& 

dV dV dV 



lim ^ ■ j dsn (s) e «»-«W *L . j dsn(s)=> \ dsh(s)=0 (39) 



The last equality follows from the fact that this term must vanish for any (non-zero) (3. The results 
in ( ) and ( ) can be seen as a simple version of the more formal statement that the boundary of a 
boundary is zero [15], i.e., the curve which bounds an area in two dimensions has no end points. 

Hopfs Umlaufsatz The results in (38) and (39) can also be related to Hopfs Umlaufsatz [6] 
which states that the tangent vector (or the normal vector or a linear combination thereof) to a 
smooth or piecewise smooth orientable curve in two dimensions rotates by 2n on making a complete 
circuit of the curve if the curve is traversed in a counterclockwise direction and by — 2n if traversed 
in a clockwise direction. For nonorientable curves, i.e., curves which crossover or intesect themselves 
it is an integer multiple of 2tt with the integer value counting the net number of signed (+/— for 
counterclockwise/clockwise) loops. This seems obvious but in fact is rather subtle to prove [6]. 

The following is an argument in favor of the Hopf Umlaufsatz but is not a proof in the strict 
mathematical sense. Consider first the piecewise smooth, actually piecewise linear, case of a polygon 
defined by the vertices v\, ■ ■ ■ ,Vn with vn + i — v\. Let l n = v n+ \ — v n . which are the directed sides 
of the polygon taken in order stepping around the polygon. We can write l n +i = A„+ii? (#n+i) 1 In 
for n = 1, 2, • • • , N + lwith A„ > and R (9) the two dimensional rotation matrix 



R(6) = 



cos [8] sin [9] 
— sin [6] cos [9] 



(40) 



Thus l n +i — A„+ii? (#n+i) • In defines each successive side vector l n as a rotated and rescaled copy 
of the previous side vector. Combining these relations for the entire polygon gives 

In+i = Atv+iAw ■ ■ ■ A2-R (#jv+i) R (On) ■ ■ R (#2) h 
= Atv+iAat • • • A2-R (9n+i + 9n + ■ ■ ■ + 8 2 ) h 

= h (41) 



which demands that 



Ajv+iAjv ' ' ' ^R (9n+i + @n + • • • + #2) — 1 (42) 



and so Ajv + iAat • • • A 2 = 1 and 8n+i + 8n + ■ ■ ■ + 8 2 = integer x (2ir) . This argument does not by 
itself specify the integer value as being unity. The fact that the integer needs to be 1 follows from 
the requirement of the polygon being orientable. 

The same type of argument can be applied to the smooth curve case. Write the tangent vector 
i at distance s along the curve as 

i (s + ds) = R{d s 8 (s) ds) ■ i (s) (43) 
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where 9 (s) is the angle between the tangent vector at s and the tangent vector at s = so that 

i(s) = R(6(s))-t(0) (44) 

Then using the fact that the curve is smooth (d s i(s) is everwhere defined) and closed we have, with 
L the total distance around the curve, that 

i(L) =R(6(L)) ■ i(0) 

= t(0) (45) 

and so 9 (L) must be an integer times 2tt. Again this result must be combined with the fact that 
the curve is orientable to specify the integer value as unity. 

3.2.2 m = 2 term 

The m = 2 term on the right equals the n — term on the left and so we have 

. , JV 1 , 



m (i) = -~ $± ■ - Z • *w) p • 

" n=l p=0 
' 1 N 



n=l 

N 



n=l 



But by definition M (1) is the area of the polygon. Defining l n = v n+ \ — v n = the vector from 
vertex n to vertex n + 1 and c„ = (v n +i + ^n) /2 = position of the center of side n we have 



Mi 



(1) = Area = ^ ($ ± ■ l n ) (p ■ c n ) (47) 

n=l 

Taking /3 = (0, 1) gives /3 ± = (-1, 0) and so 

N 

Area = - Z„,ic„, 2 (48) 

71=1 

The geometric interpretation of this is straightforward. Each of the l n ,ic n ,2 terms correspond to 
the area of a 4-sided polygon of width |Z n) i| in the x\ direction and mean height \c n $\ in the X2 
direction. If we take all the vertices to lie in the first quadrant, v n ^ > for i = 1,2 then all the 
c n ,2 are positive but the l n ,i change sign depending on whether l n points generally in the +x\ or 
— x\ direction. If we consider that the n = 1, • • • ,N ordering corresponds to following the vertices 
in a counterclockwise direction around the polygon then the positive Z nj i will generally run along 
the bottom sides of the net polygon and the negative l n ,i will generally run along the top sides of 
the net polygon, then the area of the net polygon is the total area of all the 4-sided polygons along 
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the bottom of the net polygon subtracted from the area of the 4-sided polygons along the top of 
the net polygon. 

Taking 1/2 the sum of (47) with $ = (0,1) and with $ = (1,0) we find that the area can 
also be written as 5 X)n=i det [M„] where the elements of the 2x2 matriices M n are given by 
M nt ij = v n ,iV n+ ij, which is the standard result [19]. 



Area Element We now show that the area element dv of a 2 dimensional surface dV defined by 
R(s), which is a vector in 3 dimensional Euclidean space and s = (si,s 2 ) labels points on dV, is 
given by 



dv 



dot 



d Si R(s) ■ d Sj R(s) dsids 2 = ^/ |det [ii (s) ■ ij (s)] \dsids 2 



(49) 



The two vectors dRi = d Si R(s) dsi = ii (s) dsi for i = 1, 2 (no sum on i) form a parellelogram with 
area dv whose vertices are v\ = (0,0) , v 2 = dR\, V3 = dR\ + dR 2 , v^ = dR 2 and v$ — v\ all of 
which lie on dV. (For simplicity we have assumed that Si has units of length so that d Si R(s) = U (s) 
with ti (s) automatically normalized to unity.) The area of this parallelogram dv can be computed 
directly from (46). 



dv = 



p ± ■ dR^j (p ■ dR^j + (fr ± ■ di? 2 ) (fi ■ 
j3 ± ■ -dRi) (p ■ (dRx + 2dR 2 ) 



2dRi + dR 2 



-dR- 



i? 



■dW 



= (p-e- di? 2 ) (P ■ dR^j - (p-s- <£Ri) {p ■ di? 2 ) 
= ((£■£■ ^) (P ■ ti) -(£■£■ *i) (P ■ h)) d Sl ds 2 



(50) 



The vector $ can be written as an arbitrary linear combination of the unit vectors ii (s). Letting 
$ = t\ gives 



dv = i\ ■ e ■ t 2 ds\ds 2 

= (*i, 1*2,2 - ^1,2*2,1) dsids 2 



(51) 



where tj i0 are the components of ii in an orthonormal coordinate system indexed by a = 1,2 erected 
on dV at s. But 



det [ii ■ ij] \ds1ds2 



det 



t\ ■ i\ i\ ■ t 2 
t 2 ■ t\ t 2 ■ t 2 



ds\ds 2 



\J (ti-ti) (i 2 -i 2 ) - (ii-i 2 y 

(£1,1*2,2 - £1,2*2,1) dsids 2 
dv 



dsids 2 



(52) 



For the purposes of the above derivation even though t\ ■ t\ = ti ■ ii = 1 it is more convenient leave 



dRi 



dR 2 



them as i\ ■ i\ and t 2 -t 2 . Finally note that from simple geometry we have dv 
where 9 is the angle between dR\ and dR 2 , and if we do use t\ -ii = i 2 -i 2 = 1 then v/|det [ii ■ ij] 



sin (6) 
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yl - (ii -4) = y 1 — cos (0) = sin (0) and we get the same result for dv. In general it is much 

more convenient to work with |det [£j • £,] | than other forms since the inner product ii ■ ij is 
coordinate independent. 

3.2.3 m = 3 term 

The m = 3 term on the right hand side corresponds to the n = 1 term on the left and so after 
cancelling i from both sides we have 



1 N 



71=1 



+ 



It follows from the definition of M that 



M (xi) = Area x X t 

where Xi is the "center of mass" or centroid of the polygon in the i = 1, 2 directions. 
For /3 = (1, 0), j3 ± = (0, 1) we have 

1 " 

M ( Xl ^ = T\ ( W «+l,2 - V n,2) ( V l+l,l + V l,l + Wn+l.lVl) 



(53) 



(54) 



(55) 



n=l 



and for /3 = (0, 1) , (3 ± = (-1, 0) we have 



JV 



M (ar 2 ) = K+l.l ~ U ">l) ( V l+1,2 + V l,2 + V n+ l,2V n ,2) 

Explicit evaluation of 



n=i 



dx\dxiXi 



(56) 



(57) 



with i = 1 or 2 by substituting Xi = d ■ (xfxi/2) , with no sum on i, yields the same result as in 
(55) and (56). 

3.3 Moments and Shapes of Polygons 

We begin this section by using a version of (5) modified to live in the complex plane to provide 
an alternative derivation of the result of Davis which is a generalization from triangles to polygons 
of the Motzkin-Schoenberg formula [11]. This result is also related to the so-called "shape from 
moments" problem which is to find the ordered vertices of a polygon given an appropriate set of 
the polygon moments [4] [5]. As shown by Milanfar [4], in the complex plane with z = x + iy, the 
result of Davis can be written 



J v dxdyd z 2 h(z)= l -f2 



Z n-l Z n Z n Z n+l 



_^ \ z n— 1 z n %n %n-\-l 



h(z n ) 



(58) 
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Here V is a simply connected orientable polygon with vertices z n = x n + iy n (z* = z n — iy n ), 
n = 1, • • • , N, in the complex plane and the function h (z) is analytic (= holomorphic = regular) in 
the closure of V. In the sum we have let zo — zn and z\ — zn+i- For the remainder of this section 
we follow the standard notation for variables in the complex plane, i.e., we replace xiwith x, xi 
with y, etc.. 

To write h (z) as a Fourier transform start with the definition of an analytic function, i.e., that 
it can be written as a power series in non-negative powers of z 

oo 

h( z )^j2 a ^ n ( 59 ) 

n=0 

Now define the function h ((3) , with (3 real, by 

a n = J dph(p)^- (60) 

This may seem restrictive but if we take h (/3) to vanish outside |/?| < 1 then we can represent h ((3) 
as 

oo 

h(0) = J2 A nPn(P) (61) 

n=0 

where P n (f3) are Legendre polynomials. With this representation the complex coefficients A n can 
be chosen to satisfy (35). Again assuming we can interchange sums and integrals at will we can 
write 

h (z) = [ d/3h ((3) £ { ^f- = f d(3h [fi) e*' (62) 
Using the obvious notation (a, b) ■ (x, y) = ax + by, we now have 



dxdyd z 2 h (z) = - [ d/3h (J3) [ dxdyp 2 e if}z 
v J Jv 

l - J df3h {(3) dxdy (d x , d y ) ■ (3 (1, -%) exp [i(3 (1, i) ■ (x, y)} 
If- ^ /3(1> -i)-e- ((x, u) - (x,y)„) 

x (exp (i(3 (1, i) ■ [x, y) n+1 ) - exp (i/3 (1, i) ■ (x, y)J) 



= i [ d (3h((3) y ( z * n < - 1 - 



_ exp (0z n ) 

n=l \ n n *"n+l / 

In the third line e is the Levi-Civita matrix defined in (14). In the second step we have used (5) 
and in the last step the definition of h (z) in terms of its Fourier transform. The /3's cancel in 
the coefficient of the exponents in line three and hence we reproduce the result of Davis that the 
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integral over a polygon of d 2 h (z) is the sum of h (z) evaluated at the vertices of the polygon times 
coefficients which depend only on the vertices and not on h (z) . 

The moments of a polygon obviously contain the polygon shape information. The order of the 
vertices is important as reordering them nominally leads to a polygon with a different shape and/or 
can make it nonorientable. For N vertices there are 2N independent real numbers corresponding to 
the Vi, ■ ■ ■ ,vn vertices which define the polygon. Thus the infiinite set of all possible moments of 
the polygon must be highly redunant.We will not solve this "shape from moments" problem here. 
Milanfar and others [4] [5] have shown how to solve for the vertices given a particular set of complex 
moments which can be easilly computed from (63) by letting h(z) = z fc .Here we merely present 
the complete set of relations between all possible moments M (x a , j/ b ) and the vertices (x, y) n of a 
polygon that follows from (34). 

To derive an explicit relation between the moments and the vertices of a an arbitrary orientable 
polygon multiply (32) through by j3 2 = 0\ + /3 2 and use the fact that the m = 1 term on the right 
hand side vanishes identically. The derivation is facilitated by defining a function 6 (• • • ) which 
vanishes if any one or more of its arguments is negative and equals 1 otherwise. This function can 
be used to keep track of the limits on the sums when the summation indices are redefined so that 
the powers of fi l and (3 2 are written as /3 1 a /3 2 b on both sides of the equation. Then using the fact 
that the coefficient of (3ift 2 on the left must equal the coefficient of Pil3 2 on the right for the 
same given non-negative integer values of a and b we find, with v n — (x l rl , x 2 ,„) 

M(x, a - 2 x 2 b ) , M (x-, a y 2 ~ 2 ) 



{a-2)\b\ v ' a\{b-2)\ 



= 0(6-1,0 + 6-2)^ 



a+ -- q( 0(p + q -a) (a l b 



(a+b)\ q\(a+b-l-p-q)\ (a-q)\( V +q-a)\ 
q=0 p=0 \ x Z^n=l l x l,"+l x l,"Ml,n+l x 2,n+l x l,n x 2 
la+b-l-q / g(„j_„j_i „\ 1 (a+6-l-p)! p 



-6 (a -la + &-2lV V \ e (P + 1 + 1 a ) ^+&yi q \( a +b-i-p-q)\ ( -g-i)!(p+ 9 +i-o)! 

U(CL 1,0-1-0 £) 2^ xV W fr, x,-T, W« a+b-l-p-q a-q-1 p+q+l-a 

(64) 

The derivation is straightforward but tedious. 



3.4 Fraunhofer Diffraction 

Fraunhofcr diffraction occurs when the light diffracted from an aperture or opening in an opaque 
screens is observed in a plane far from aperture itself [2]. By "far" we mean that the distance 
between the opaque screen and the plane of observation, which is by convention taken to be parallel 
to the screen, is much larger than the maximum dimension of the aperture. For the case where 
the opaque screen lies in the X\x 2 plane and the illumination is a unit amplitude plane wave of 
wavelength A incident on the screen from the side opposite the observation plane, the amplitude 
of the diffracted light at position a;, in the observation plane a distance L away from the screen is 
given by 

A (a?) = J d 2 x'6 v (x 1 ) exp 

v 



1—XiX; 

L % 



d 2 x' exp 



1—XiX: 

L 1 



(65) 
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Here Qy (x,y) = 1 inside the aperture and outside describes the aperture shape and k = 2tt/X. 
This is essentially the left hand side of (5) but with 



_ kxi 
1 l ~~L 



(66) 



and without the normalization factor. The intensity of the diffraction pattern is given by / (x) = 
\A(x)\ 2 . 

We consider two cases, a circular aperture whose solution follows almost trivially from (5) and 
a slit whose solution is already implicit in (32). 

For a circular aperture of radius R centered at Xi — it follows from (5) that 



y ' ip 1 



2tt 



dipr (ip) exp i/3 • Rf (ip) 



(67) 



with f (ip) — x'/ \x' | . Writing j3 ■ f (<p) — (3 cos {ip — ip^) where <pp is the angle (3 makes with respect 



to the xi-axis and 



P 



ksjx 1 + y 2 /L, (67) becomes 



A(x) = 



R 

ip 



2tt 



dip cos {ip — ipp) exp [i(3Rcos (ip — <Pp)\ 



271 



R — i f 
= -^-rd R / dip exp [i (3 R cos (<p)} 
ip P J 
o 

= -^2wd R J (pR) 

- '* ^) 



(68) 



Here the J n (x) are the Bessel functions of the first kind. This is the standard result for a circular 
aperture [2]. . 

For a slit (rectangular) aperture of width 2a^ in the Xi direction, centered at Xi — 0, we merely 
have to substitute v\ — (a\, ai) , v 2 — (— Q-i, 0,2) , «3 = (— ai, —a 2 ) , U4 = (01, —02) , and V5 = v\ into 
(32) to obtain 

A(x) = 



N 



>_ P ■ £ ■ {Vn+1 ~ Vn) 

P 2 n=1 P ■ (v n +i - V n ) 



(exp (iP ■ Vn+i'j - exp (ifi ■ v n ^j ) 



P 2 



( (exp (—iaiP 1 + ia 2 P 2 ) - exp {iP x ai + iP 2 a 2 )) \ 
(exp (-ia 1 p 1 - ia 2 P 2 ) - exp (-ifi^i + iP 2 a 2 )) 
+^ (exp (id! Pi - ia 2 p 2 ) - exp (-iP^t - iP 2 a 2 )) 
y (exp (ioi/3! + ia 2 p 2 ) - exp(+iP 1 a 1 - iP 2 a 2 )) J 

sin ( / 5 1 ai) sin (P 2 a 2 ) 



(2ai) (2o 2 ) 



(69) 
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which again is the standard result [2]. 

Fraunhofer diffraction patterns for arbitrary (orientable) polygons can be calculated simply by 
substituting the vertex values into (32). It is interesting to compare the patterns generated for a 
given set of vertices as the vertices are reordered to make the polygon nonorientable. 



4 Results in Three Dimensions 

Before discussing Porods law we point out an essentially obvious but useful fact. Since the faces 
of a polyhedron are themselves polygons it follows that applying (5) to a polyhedron reduces the 
integral over the volume to a sum of integrals over the areas of polygonal faces and then applying 
(32) to the faces themselves reduces the volume integral to a sum of integrals over the edges. These 
integrals can of course be evaluated exactly. The only issue is the bookkeeping required to keep 
proper track of the vertices. 



4.1 Volume Element 

We use (5) to show that the volume in three dimensions of a parallelpiped formed by three vectors 
di is given by det [ay] where ay the 3x3 matrix of the components of the a*j, i.e,.the first column 
is an, the second an and the third a 3 i with of course, i = 1,2, 3. If the a*j are infinitesimals given 
by dRi = d Si Rdsi with no sum on i, then det [dRij] d 3 s forms the fundamental volume element dv 
for integration. Expanding the right hand side of (5) to first order in f3 we have that the volume of 
a shape in three dimension is given by 



</■•■ 

v 

= ft- J d 2 s^fi(i)n{s)(p-R{s)) (70) 

OV 

A parallelepiped, defined by three (non coplanar) vectors a*j, i — 1,2,3, in three dimensions has 6 
parallelogram faces. The functions Rf (s) for positions on the faces with / labeling the faces, are 

R\ (s) = si&i + s 2 d 2 R 2 (s) = sid 2 + s 2 s 3 R 3 (s) = Si&s + s 2 ai 

Ri(s) = d 3 + sidi + s 2 d 2 R 5 (s) = a*i + sid 2 + s 2 s 3 R e (s) = a 2 + sia 3 + s 2 di (71) 

The values of Si in each case range from to the corresponding length of the associated vector 
di. For example, for / = 3, si ranges from to 1 03 1 = a 3 and s 2 ranges from to |a*i| = 01. We 
assume that the vectors di are ordered so that (5*2 x 0*3) points in the general direction of a*i, i.e., 
di ■ (a 2 x 03) > where "x" is the standard cross product. On each face <?/ (s) and n/ (s) are 
constants and so factor out of the integrals. The outward normals fif on the opposite faces, 1 and 
4, 2 and 5, 3 and 6, point in opposite directions. Thus, the integrals over the s dependent parts of 
Rf (s) completely cancel and we are left with 



v = 



(ft ■ ^/gln^j (j3 ■ a 3 ) aia 2 
+ (P ■ V95n 5 ^j (jl ■ a^j a 2 a 3 

+ (ft ■ v^^e) (ft ■ a 2 ) a 3 «i (72) 



17 



But $ is arbitrary. Choose it to be &i. Then since fi\ = (d\ x a 2 ) / |fii x a 2 | and n 6 = (S3 x fii) / |&3 x Si| 
we have a\ ■ h 4 = a\ ■ n 6 — and u reduces to 



v = (ai ■ y/gEn 5 ) axa 2 a 3 



(73) 



But, as shown above in the derivation of the area element in two dimensions, y/g is the sine of the 
angle between the two correponding vectors and so equals the magnitude of the cross product of 
those vectors, hence y/gE = \a 2 x a 3 | and thuSy/g^n 5 = a 2 x a 3 .So finally 

v = d\ ■ (<Z2 x as) a 1 a 2 a 3 
= «i • a 2 x i 3 

= £ijka>liCl2jCl3k 

= det [a,-,-] 



Met [ai ■ dj] (74) 
Again the form in the last line is the most useful since it is coordinate independent. 

4.2 Porods Law 

Finally, we rederive the anisotropic version of Porods law as given in the work of Ciccariello, et. 
al. [12]. The anisotropic result of course reduces to the isotropic result for spherical particles. Here 
we use the notation k rather than ft as is more common in this context. 

The intensity / (jtj — I {kk^j = I (x^) °^ light of wavelength A scattered off a particle defined 

by the shape function 9y (x) in the direction k = k out — ki n where ki n (jt ou tj is the direction of 
propagation of incident (scattered) light is 



CO 



d D x exp ik ■ x 



(x) 



(75) 



To be general, for now, we start in D dimensions. Porods law [12] follows from this in the case 
where the magnitude of the scattering wavevector k = k = 2ir/\ is large, i.e, A is small and so 
Porods law is often associated with X-ray scattering. 

To obtain Porods law, use (5) and evaluate the integral over dV in the limit of large k using the 
method of stationary phase. The positions a p on the surface R (s) where the phase is stationary 
satisfy 

kd s (k-R(s)) =0 (76) 



d, 



(k-R(s)) 



Here d s — (d Sl , d S2 , • • • ) is the gradient with respect to the surface coordinates and the index p runs 
from 1 up to the number of solutions of (76). 

For simplicity consider a single solution a to this equation. The argument of the exponential in 
the integrand can now be approximated by 



ik ■ R(s) = ik ^fc • R{a) + i {d ai d aj k ■ R(a)j (s, - <Ji) (sj - 



(77) 
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to second order in s — a. Substituting into (5) gives 



d D x exp 



v 

~ exp 
~ exp 



ik ■ , 



ikk ■ R (a) 
ikk ■ R (<?) 



k_ 

ik 



exp 



ik 
~2 



(dendcTjk ■ R(v)^J (si - o-j) (sj 



ik 



h(a) 



det 



-%d ai d aj k-R(a) 



(78) 



where the determinant in the denominator is taken over the i,j indices. Using the fact that 



det 



- l M ai d aj k-R(a) 



= fc^-Met 
d D x exp 



-id ai d aj k-R(a) 



we get 



X 



and so 



ik ■ x 



I[k 



1 



1 



1 

fcfc^- 1 )/ 2 ~ fc(^+i)/ 2 
1 



k D+l 



(79) 



(80) 



Thus in 3 dimensions / ysj scales as 1/fc 4 for large k. This is the standard isotropic statement 

of Porods law [3]. But as pointed out by Ciccariello, det d ai d aj k ■ R(a) is proportional to the 

Gaussian curvature (in 3 dimensions) of the surface dV. In general for an arbitrary shaped particle 
there will be multiple stationary phase points, i.e., <?i,<72, • • • which must be summed to get the 
complete amplitude. There are also issues with positive and negative curvature which must be 
carefully considered as well along with how the handle the case when the curvature vanishes [12]. 



5 Conclusion 

We have shown how a simple idea, that of combining Gauss's law with the Fourier transform 
provides alternative solutions and/or derivations of many different classical results in physics and 
mathematics. No doubt there are many other problems and proofs to which this idea can be applied. 
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